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Abstract 

Approximate probabilities of formation of various atoms and ions in different bound states are 
determined for the exothermic nuclear (n, ^Li; t, ^He)-reaction of atomic lithium-6 with slow neu- 
trons. In our calculations of the final state probabilities we have assumed that the incident lithium 
^ ■ atom is in its ground (doublet) atomic '^S{L = 0)— state. It is straightforward to generalize our 

^ ^ \ analysis to other bound states of the three-electron Li atom. 

> 

oo 
en 

o 



% 



* E-mail address: lafrolov@uwo . cal 
'''E-mail address: ldwardlaw@uwo . cal 



I. INTRODUCTION 

The nuclear reaction of ^Li nuclei with slow neutrons is written in the form [l| 

^Li + n = ^He + t + 4.785 MeV (1) 

where the notations ^He and t stand for the helium nucleus (often also called the a— particle, 
or a, for short) and tritium (or ^H) nucleus. For thermal neutrons with En ^ Q the cross- 
section 0" of this nuclear reaction is very large (Jmax ~ 960 ■ 10^^^ cm? or 960 ham [2], for 
short. The velocities of the two nuclear fragments formed in the reaction, Eq.([T]), with slow 
neutrons are Vt ~ 6.03986 a.u. and Va ~ 4.52989 a.u. for the tritium nucleus and a— particle, 
respectively. These velocities are given in atomic units, where % = l,me = l,e = 1. The 
unit of atomic velocity is Ve = ac ^ ^ ^ 2.1882661 ■ 10^ cm ■ sec~^, where c is the speed of 
light and a is the fine structure constant. Formally, the atomic velocity Ve is the velocity of 
the Is— electron in the hydrogen atom with the infinitely heavy nucleus °°H. It is clear that 
in atomic units Vg = 1. 

The nuclear reaction of the ^Li nucleus, Eq.([T]), proceeds with neutrons of all energies and 
the energy released increases almost linearly with the energy of the incident neutron. For 
fast neutrons with £"„ > 1 MeV the energy released is substantially different from the value 
quoted in Eq.(Il]) and the velocities of the atomic fragments increase correspondingly. Note 
that the cross-section of this reaction has a large resonance (maximum) a f^ 4.5 barn at 
En ~ 240 — 270 keV, but it is relatively large for neutrons of all energies En < 0.8 MeV. This 
makes the reaction, Eq.(IT]), extremely important in the thermonuclear ignition and following 
propagation of the thermonuclear burning wave in highly compressed (p > 100 g ■ cm,~^) 



^LiD deuteride which is routine 



used as a thermonuclear fuel in modern thermonuclear 



explosive devices (see, e.g., |3[, ^). The analogous (n, ^He)-reaction [S] also plays a very 
important role in such processes. 

In general, the nuclear (n, t)-reactions of the ^Li and ^He nuclei with neutrons allow one 
to: (1) reduce drastically the overall bremsstrahlung loss from the combustion zone; and (2) 
increase the tritium/deuterium ratio which is crucially important to start many new ther- 
monuclear ((i,t)— reactions. Briefly, by using the ^LiD deuteride in modern thermonuclear 



explosive devices we can reduce the required compressions to relatively small values. 



in 



many cases such compressions are dozens times smaller (usually, in 25 - 40 times smaller j^) 
than compressions required for any other (solid) thermonuclear fuel, e.g., ^LiD deuteride. 



On the other hand, by using compressions which are provided by a standard 'primary' nu- 
clear charge one can create extremely compact thermonuclear explosive devices based on the 
^LiD deuteride. The idea to use the pure ^LiD deuteride in thermonuclear explosive devices 
was originally proposed by V.L. Ginzburg in 1948-1949 (see discussion and references in |3|). 

Our goal in this study is to make accurate numerical predictions of the final state prob- 
abilities for the reaction, Eq.([T]), with slow neutrons. We want to predict (accurately) the 
probabilities of formation of various atoms and ions in different (ground and excited) final 
states. Briefly, we want to determine the probabilities to detect the He and ^H atoms and 
He"*", ^H~ ions in their bound states. Note that all newly created atomic fragments from re- 
action, Eq.([T]), move rapidly even in the case of slow neutrons and that accurate theoretical 
prediction of the final state probabilities for rapidly moving atomic fragments is not trivial. 

On the other hand, the results of our evaluations and methods created for such evaluations 
are of great interest in many applications related with the nuclear reaction, Eq.([T]), and 
another similar reaction of ^°B nuclei with slow neutrons, namely 
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B + n = ^Li + ^He + 2.791 MeV (2) 



The reaction, Eq.(|2]), is extensively used in the 



Doron neutron capture therapy (BNCT, for 



i-Q 



to treat different forms of cancer. 



short), or boron neutron-capture synovectomy 
including brain cancer. The fast a— particle produced in the reaction, Eq.([2]), kills (or at 
least 'badly damages') one cancer cell before it finally stops. The modern applications of this 
reaction to cancer treatment are based on the use of molecules which contain a large number 
of ^°B-atoms, e.g., the NasfBon HivNHo CHoCHo NHol molecule, other similar molecules, and 

n n n ^ ^ ^^ ' 

molecular clusters [8|, [9[ (see also p^ and references therein). In this case the overall energy 
release from the reaction, Eq.([2]), in one cancer cell can be extremely large. Correspondingly, 
the local temperature in the whole cell suddenly increases to very large values and this kills 
the incident cancer cell with almost 100 % probability. 

Note that the tritium nucleus does not form in the nuclear reaction Eq.([2]) which means 
that it is safe to initiate this reaction inside of a human body. By studying the nuclear 
reaction, Eq.(Il]), in few-electron atoms and ions we want to develop a number of reliable 
theoretical methods and numerical procedures which can be later used in applications to the 
analogous reaction, Eq.([2]). 



II. APPROXIMATE VARIATIONAL WAVE FUNCTION OF THE LITHIUM 
ATOM 



In our earlier paper [5| we have calculated the final state probabilities to form one-electron 
atomic species in those cases when exothermic nuclear (n;t)— and (n; a)— reactions occur 
in one electron atoms and ions. In this study we deal with the actual three-electron wave 
function of the Li atom and two-electron wave function of the He atom. For our purposes in 
this study it is important to construct accurate variational wave function(s) of the ground 
(doublet) ^S{L = 0)— state of the Li atom which is written in the following general form 



(see, e.g., (ll|, (l2|) 



\I^(Li)j;^=o = ipL=oi^] {rij}){af3a — /3aa) + (pL=o{B; {rij}){2aaj3 — j3aa — a/3a) (3) 

where 'iIjl=q{A] {ry}) and (J)l=q{B] {vij}) are the two independent radial parts (= spatial 
parts) of the total wave function. Everywhere below in this study, we shall assume that 
all mentioned wave functions have unit norm. The notations a and j5 in Eq. ([3]) are the 



one-electron spin- up and spin-down functions, respectively (see, e.g., ISj). The notations 
A and B in Eq.dH]) mean that the two sets of non-linear parameters associated with radial 
functions -0 and can be optimized independently. Note that each of the radial basis 
functions in Eq.([3]) explicitly depends upon all six interparticle (or relative) coordinates 
^"12, '"is, '"23, ru, r24, r34, where the indexes 1, 2, 3 stand for the three electrons, while index 4 

mean the nucleus. 

In our earlier work [M] we have introduced an advanced set of radial basis functions 
for bound state computations of three-electron atomic systems. Such a set is called the 
semi-exponential basis set and it is written in the form 

N 

fc=i 
where a^, /3fc, 7^ {k = 1,2,..., A^) are the varied non-linear parameters. The use of a large 

number of non-linear parameters in Eq. (jl]) allows one to construct compact and very accurate 



variational wave functions for different three-electron atoms and ions. It was shown in [11 1 
that the semi-exponential basis, Eq.(j3]), has a large number of other advantages in accurate 
numerical computations. 

In the sudden approximation [l5|, [l6| the numerical determination of the final state 
probabilities is reduced to the analytical computation of the Fourier transform of the overlap 



integral between the incident ipin and final ifjfi wave functions. In the case of the nuclear 
reaction of ^Li nuclei with slow neutrons under consideration here, the incident wave function 
represents the ground ^S*— state of the lithium atom. The final wave function represents the 
product atom and/or ion. For the helium atom product we need to determine the following 
integral (or probability amplitude Mjj) 



Mif = ^Li(ri4, r24, r34, ri2, n^, r23) exp(zV„ ■ r^ + zV^ ■ r24) x (5) 

^He(?^14, r24, ri2)dhudh24dh34 

where \E'Li is the wave function of the Li atom, while \t'He is the wave function of the final 
(bound) state of the He atom. The velocity V^ is the final velocity of the a— particle 
after reaction, Eq.(IT]). In reality, we need to determine the correspponding probabilities 
Pif = M*p-Mif. This allows one to reduce the problem (see below) to numerical computation 
of the one- and two-electron density matrix. 

For the tritium atom product in one of its bound states the probability amplitude Mif is 
written in a slightly different form 



M 



if 



J J J ^Li(ri4, r24, r34, ru, na, r23) exp(zVt ■ ri4)^H(ri4)<i^ri4(i^r24(i^r34 (6) 



where \&h(^14) is the unit-norm wave function of the tritium (or hydrogen) atom and Vt is 



its velocity after the reaction, Eq.(IT]). In our earlier study [17(| we have calculated integrals 
similar to the integrals in Eq.(|5]) and Eq.([6]) for the two-electron atomic systems which are 
involved in the analogous nuclear reaction of the ^He nuclei with slow neutrons. The integral, 
Eq.(|5]), can be considered as a partial Fourier transformation of the overlap of the incident 
and final wave functions. This integral represents the Galilean transformation between the 
incident system (which was at rest) and the final system which is rapidly moving with 
constant speed. 

In general, numerical calculations of integrals Eqs.(|5]) and (|6]) are difficult to perform due 
to the presence of the electron-electron coordinates Vij in the wave function. On the other 
hand, for approximate evaluation of the final state probabilities for the nuclear reaction, 
Eq.(Il]), we do not need to use the highly accurate wave functions of the Li and He atoms. In 
particular, we can use the wave functions of few-electron atomic systems which do not contain 
any of the electron-electron coordinates. This drastically simplifies analytical/numerical 
calculation of the integrals, Eq.([5]) - Eq.([6]). Formally, such approximate wave functions 



depend on the electron-nucleus coordinates only. Therefore, these wave functions correspond 
to a system of A— electrons which do not interact with each other and can be considered 
as free-electron wave functions. For the He atom we have A = 2, while for the Li atom we 
have A = 3. The explicit construction of the approximate free-electron wave functions for 
few-electron atomic systems is discussed below. 

A. Special form of the trial wave function of the lithium atom 

To avoid problems related to the analytical computation of the Fourier transformations 
given by Eq.(|5]) one can apply approximate variational expansions of the three-electron wave 
function. To simplify all following calculations we represent the wave function of the Li atom 
in the form 

ipL=o{riA, ^24, r34, 0, 0,0) = ^ Cfcr^'^''V^'^''V^'^''^exp(-afcri4 - /3fcr24 - 7fc^34) (7) 

fc=i 

Ns 
k=l 

where Ck are the linear (or variational) coefficients, while mi{k),m2{k) and rrislk) are the 
three integer (non-negative) parameters, which are, in fact, the powers of the three electron- 
nucleus coordinates rj4 = Tj (i = 1, 2, 3). Below, we shall assume that the trial wave function 
Eq.([7]) has a unit norm. Furthermore, in all calculations performed for this study only one 
spin function xi = aj3a — j3aa is used. It is clear that the wave function Eq.([7]) contains only 
electron-nuclear coordinates and does not include any of the electron-electron coordinates. 
The real (and non- negative) parameters ak-,Pki1k are the SA^^ varied parameters of the 
variational expansion, Eq.([7]). The wave function, Eq.([7]), must be properly symmetrized 
upon all three electron coordinates. This problem is discussed in the next Subsection. 

The principal question for the wave function, Eq.([7]), is related to its overall accuracy. If 
(and only if) such accuracy is relatively high, then such a wave function, Eq.(I7]), can be used 
in actual computations of the probability amplitudes, Eqs.(|5]) and ([6]). In actual applications 
the approximate wave function, Eqs.(j5]) and (J6]), can be constructed from the highly accurate 
wave functions already known from earlier works (see, e.g., |ll( and references therein). 
Briefly, we can take our wave function of the ground ^S"— state of the Li atom from [llj] and 
remove all those terms which contain electron-electron r^j coordinates. Then the non-linear 



parameters in the trial wave function, Eq.(I7]), must be re-optimized. The resulting wave 
function can be considered as an optimal independent-electron wave function of the ground 
state of the °°Li atom. Using this approach we have determined the 23-term variational 
wave function shown in Table I. The total energy E of the ground ^S*— state of the °°Li 
atom obtained with this independent-electron wave function is -7.44859276608 a.u. Note 
that this value of E is close to the exact total energy of the ground state of the °°Li atom. 
This indicates a very good overall quality for our approximate wave function which does 
not include any of the electron-electron coordinates ri2,ri3,r23. Note also that in atomic 
physics based on the 'Hatree-Fock' and even 'hydrogenic' approximations the ground state 
in the Li atom is designated as the 2^5*— state, while in the classification scheme developed 
in highly accurate computations the same state is often designated as the l^S"— state. This 
classification scheme is very convinient to work with trully correlated few-electron wave 
functions. It is clear that no hydrogenic quantum numbers are good in such cases, and we 
have to use the more appropriate (and convinient) classification scheme. To avoid confllicts 
between these two classification schemes in this study we follow the system of notation used 



by Larsson [12| which designated this state in the Li atom as the 'ground ^S'-state'. 



B. Antisymmetrization of the trial wave function of the lithium atom 

The actual many-electron wave function in an atomic system must be completely an- 
tisymmetric with respect to all electron variables, i.e. upon all electron spatial and spin 
variables. This statement is true for all exact and approximate few-electron wave func- 
tions, including the optimized free-electron wave functions, Eq.([7]). Antisymmetrization of 
the two-electron wave function is trivial and is not discussed here. For a three-electron 
atomic wave function this requirement is written in the form ^123^1/(1,2,3) = — \1/(1,2,3), 
where \1/ is given by Eq.(l3]) and Ae is the three-particle (= electron) antisymmetrizer 
Ae = e. — Pi2 — Pi3 — P23 + -Pl23 + Pi32- Here e is the identity permutation, while Pij 
is the permutation of the i-th and j-th particles. Analogously, the operator Pijk is the per- 
mutation of the z-th, j-th and fc-th particles. In actual computations antisymmetrization of 
the total wave function is reduced to the proper antisymmetrization of corresponding matrix 
elements (for more detail, see, e.g., [ij]). Each of these matrix elements is written in the 
form (^ I O I \1/), where O is an arbitrary spin-independent quantum operator which is 



symmetric upon all interparticle permutations. The wave function \1/, Eq.([3]), contains the 
two different radial parts i(j and (p. By performing the integration over all spin coordinates 
one finds the four following spatial projectors V^^,V^^ = V^^p and V^f,^ 

V^^ = ^ (2e + 2 A2 - A3 - A3 - A23 - A32) (8) 

V^^, = - i^Pis — P23 + -P123 — -P132 j (9) 

n^ = I (A3 - A3 + A23 - A32) (10) 

V^^ = ^ (2e - 2A2 + A3 + A3 - A23 - A32) (11) 

Here the indexes "0 and (p correspond to the notations used in Eq.(j3]) to designate the 
two spatial parts of the total wave function. For an arbitrary symmetric spin-independent 
operator O each of these four projectors produces matrix elements (\1/ | O | \1/) of the correct 
permutation symmetry (for doublet states) between all three electrons. 



C. Bound state wave functions of the final atomic fragments 

The final atomic states arising in the exothermic nuclear (n, ^Li; t, ■*He)-reaction of atomic 
lithium-6 with slow neutrons contain either one, or two, or zero bound electrons. In this 
study our main interest is restricted to one- and two-electron atoms and ions. The explicit 
form of one-electron atomic wave functions takes the form (see, e.g., [l8|) $„£m.(^) ©) <P)<^ = 
RneiQ, r)Yim{Q, 4>)a, where a is the spin-up wave function, ^£^(6, 0) = Yfm(n) is a spherical 
harmonic and Rne{Q, r) is the radial function. The radial function is written in the form 



1 Q(n-£-l)!r2Qrnm"^-i(-l)M n + e \ .2Qr.k . Qry 



rn 




n 



n 



where Q is the nuclear charge, while n and I are the quantum numbers of this bound state. 
Note that the radial function, Eq. (TT2|) . has a unit norm for an arbitrary Q. 

The wave function of the two-electron He-atom can be approximated in a variety of 
different forms. Currently, it is possible to construct approximate wave functions which 
provide 25 - 50 correct decimal digits for the total non-relativistic energy of the He atom 
and He-like ions. In this study, however, we shall use an approximate He wave function 
which does not contain the electron-electron coordinate r2i- For the ground l^S"— state of 
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the °°B.e atom the radial part of such a wave function is written in the form 

r 5 5 1 

^L=o(ri,r2,ri2) = Cexp -{Q - —)ri - {Q - —)r2 (13) 

L io iO J 

(Q— — ) 
where Q = 2 is the nuclear charge for the He nucleus and C = — ^^^ is the normal- 
ization constant. The corresponding spin part of the total wave function takes the form 
rj = a{l)l3{2) — /3(l)a(2) = a/3 — /3a. The explicit form of the rj spin function is important 
in the performance of integration over all spin coordinates. 

The total energy of the °°He atom obtained with this wave function, Eq. flT^ . is -2.85 
a.u. which indicates that it provides a relatively good approximation to the actual ground 
state wave function. The main advantage of the approximate wave function, Eq. (TT3|) . is 
its explicit dependence upon the two electron-nuclear coordinates ri and r2 only. This 
drastically simplifies the subsequent computation of radial integrals with the Bessel functions 
(see below). On the other hand, we need to note that the best-to-date one-term radial wave 
function for the °°He atom (see, e.g., [5|) corresponds to substantially better numerical 
accuracy, since it provides the total energy E = -2.899 534 375 443 69 a.u. which is very 
close to the exact answer (see, e.g., [l9(]). The corresponding non-linear parameters can be 
found in [5|]. But, in contrast with the wave function, Eq. (ll3l) . the wave function from 
explicitly depends upon the electron-electron coordinate r2i. 



III. CALCULATIONS 

By using the free-electron wave function of the Li atom obtained above we can estimate 
the probabilities to form various atomic species during the nuclear reaction, Eq.([T]). The 
structure of the trial wave function, Eq.(I7]), enables one to perform accurate computations 
of all integrals which include one, two and even three Galilean exponents for electrons. In 
reality all such integrals are reduced to the products of one-dimensional integrals. In other 
words, in our approach all electron coordinates are separated and this simplifies drastically 
the following analytical and numerical computations of all required integrals. 

To determine the matrix element of the operator exp(2V-r) (Galilean exponent) between 
the wave functions of the initial and final bound states, we apply the Rayleigh expansion of 
a plane wave (with the use of spherical harmonic addition theorem [18|], |20|]). An alternative 
approach is based on direct calculations of the integrals containing the exp(— Ar + tV ■ r) 



factor. The explicit formula takes the form (see, e.g., 18|], 20|, 21| ) 



e=o m=-e 



(14) 



where ^^^(n) are the spherical harmonics, V is the velocity of the final atomic fragment 
and rij^ = - is the unit norm vector which corresponds to an arbitrary non-zero vector y. 
Also in this equation the spherical Bessel functions je{Vr) are defined by the relation (see, 

e.g., y , m 



MVr) = sj' 



— J 

2Vr '-• 2 



(15) 



where J^_,_i(x) are the Bessel functions. These formulas are used in analytical and/or nu- 
merical computation of all required matrix elements. 

Actual computations of matrix elements with the 'factorized' trial wave functions, Eq.(I7]) 
are performed with the use of the following formula 

T{ii + v + l) ^h\u 1 



t^Ju{bt)exp{-pt)dt 



r(zy+l) V2 



2^i( ,^^;z/ + l; 



62 



(16) 



2 ' 2 ' ■ ">2 + 52y 
where the notation 2Fi{a,b;c]x) stands for the hypergeometric function. In many actual 
cases these hypergeometric functions are related to the elementary /rational functions, since, 
e.g., 2F'i{a,b;b]x) = (1 — x)~". Note also that all calculations for this study have been 
performed with the use of one- and two-electron density matrices 2J| , 25|] constructed for 
the incident atomic systems. The one-electron density matrix for three-electron Li atom is 
defined as follows 



p(x3,y3) = J J \ ^(xi,X2,X3))(^(xi,X2,y3) I d^Xid^X2 



(17) 



The probability to find the final one-electron atom/ion in the bound state with the wave 
function (pai'^s) is now written as the following double integral 

P = J J exp(zV ■ r3)0a(x3)p(x3, y3) exp(-«V ■ r'^)(j)aiy3)d^^3d^y3 (18) 

where Xj = (rj, Sj) and y^ = (r^, s'j) are the spin-spatial coordinates of the i-th electron. It 
is assumed in Eq. OlSp that the bound state wave function (pa is real. If the one-electron 
density matrix p(x3,y3) can be represented in a factorized form, e.g., in the form 



N N 



p(x3, ys) = E E C^Cjr''^ exp(-7ir3)(r7"^ exp(-7jr^) 



(19) 
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then the formula for the final state probability is reduced to the form 

P =1 M^^f 1 2= MUfMi^f (20) 

where Mi^f is the probability amplitude. To compute this probability amplitude one needs 
to use only one-dimensional integrals. It follows from the definition of the one-electron 
density matrix that 

//p(x3,y3)c^'x3dV3 = l (21) 

In actual calculations this condition is an important test of the constructed density matrix. 
Other tests were based on computation of different bound state properties of the Li atom 
with the use of one-electron density matrix. The two-electron density matrix is defined 
absolutely analogously and we do not want to discuss it here. The two-electron density 
matrix for the Li atom is needed in those cases when the goal is to evaluate the probability 
of formation of the He atom in different bound states. 

The results of our computations of the final state probabilities for one-electron atomic 
species formed in the reaction Eq.([T]) can be found in Table II. As mentioned above these 
final states include various ground and excited states in the helium atom (^He), tritium 
atom ('^H) and He-like one-electron ion (He"^). In this study we restrict ourselves to the 
consideration of the Is— ,2s— and 2p— states in the final one-electron ^H atom and "^He^ 
ion (see Table II). For the neutral ^He atom we have evaluated the final state probability 
only for its ground l^S"- state (see Eq. (ll3p ). We have found that the corresponding final 
state probability is ~ 7.45719-10"'^ %, i.e. it is a relatively small value. Note that each 
nuclear fragment from reaction, Eq.([T]), has a very large velocity. Therefore, it is very hard 
to observe the final atomic species with bound electrons. Based on these arguments one can 
predict that the main atomic products from reaction, Eq.([T]), will be the tritium ion ^H+ 
and helium ion ^He^"*". The overall probability to form other atomic species is less than 1 %. 

IV. CONCLUSION 

We have considered the nuclear reaction, Eq.([T]), involving the ground ^S"- state of the 
three-electron Li atom. The probabilities of formation of different atomic species during this 
reaction have been evaluated numerically, by using the actual three- and two-electron wave 
functions for the Li and He atoms, respectively. Our newly developed procedure is based on 
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the use of the optimized free-electron wave functions for few-electron atomic systems involved 
in the process. This allows us to perform all required complete and/or partial Fourier 
transformations of the wave functions. The computed final state probabilities are close to 
the exact values determined with the use of highly accurate (or completely correlated) wave 
functions which include all electron-electron (or correlated) coordinates. The method used 
in our procedure is based on the use of one- and two-electron density matrices constructed 
for the ground state of the three-electron Li atom (for more detail, see e.g., j25|]). 

Our procedure can now be used for the more complicated nuclear reaction Eq.(|2]) in the 
five-electron B-atom. The Li-atom/ion and He-atom/ion which form during this reaction 
may contain up to three and two electrons, respectively. These atomic fragments move 
rapidly, with the velocities v^ ~ 2.40896 cm ■ sec~^ and Va ~ 4.21568 cm ■ sec~^ in the case 
of slow neutrons. The sudden approximation can thus certainly be applied to the He-atom 
and He-like ions. However, this approximation cannot be used for internal electrons (or 
l^s— electrons) of the Li atom/ion, since the velocities of these two electrons are comparable 
with the final velocity of the ^Li nucleus. It is very likely that the probability to observe one 
of the ^Li^"*" and ^Li"^ ions after reaction Eq.(l2]) will be relatively large (> 20 %). Note also 
that, if the nuclear reaction Eq.([T]) is produced by a fast neutron with E^ > 1 MeV, then 
the velocity of the Li atom/ion(s) is larger than the atomic velocities of the l^s— electrons 
in the Li atom and the sudden approximation can be applied. 
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TABLE I: An example of the trial, three-electron wave function constructed with the use of A^ = 23 
semi-exponential radial basis functions, Eq.(7). This wave function produces the total energy E 
= -7.44859276608 a.u. for the ground ^S*— state of the °°Li atom. Only one electron spin-function 
Xi = a/3a — /3aa was used in these calculations. 

k mi(k) m2(k) ms{k) Ck ak i^k Ik 

0.146131429481911E+02 0.416423958308045E+01 0.390330393520923E+01 0.754647967615054E+00 
0.110681883475133E+03 0.526120133598745E+01 0.699723209902291E+01 0.604874234309236E+00 
0.234004890152000E+03 0.598950989816454E+01 0.378246213278597E+01 0.590931685390717E+00 
0.116882861303726E+03 0.413820964104499E+01 0.322650638788059E+01 0.590283203974117E+00 
-0.112533836667930E+02 0.404234943769060E+01 0.402575246908481E+01 0.136870242306376E+01 
0.124148253762594E+01 0.772078877480103E+01 0.208089611437233E+01 0.906692873308709E+00 
0.129950444032051E+01 0.125634586976915E+02 0.167180833975207E+02 0.111484879880926E+01 
0.121940082563610E+02 0.336154087416427E+01 0.329089537444583E+01 0.569461303170668E+00 
-0.815264312642982E+01 0.315618697991587E+01 0.293590042487235E+01 0.642752252837086E+00 
0.231307852206996E-01 0.216119418842614E+01 0.354528921898964E+01 0.848100637437686E+00 
-0.290323204459308E+02 0.3031 11504960694E+01 0.340407714984374E+01 0.439831339492880E+00 
-0.782805087195444E+01 0.319874342329211E+01 0.294880236447181E+01 0.670201325893484E+00 
0.466451849833029E+05 0.257191421047538E+02 0.198922117452300E+02 0.745721864825035E+00 
0.102312737866967E+03 0.610735801244370E+01 0.126738573321754E+02 0.583192334409026E+00 
0.496831723313492E+01 0.214720891197913E+01 0.139151014500059E+02 0.332004403600092E+00 
0.331570408306929E+01 0.313893628304852E+01 0.288670468153774E+01 0.656255394173358E+00 
-0.224488896165048E+03 0.111494404746535E+02 0.712253949382151E+01 0.412242317869619E+00 
-0.200426722868461E+01 0.282336074847368E+01 0.292360100848728E+01 0.638644442755579E+00 
0.186252350189224E-02 0.103900277658347E+02 0.241603283627644E+01 0.858660133717424E+00 
0.314011900712871E+01 0.314575640023816E+01 0.238785024725699E+01 0.645792868876097E+00 
0.164525876141841E+01 0.318702139403640E+01 0.348200049008609E+01 0.663487519449685E+00 
0.549306690646211E+01 0.313423005590041E+01 0.376361954006700E+01 0.494738425712773E+00 
0.489167893642778E+01 0.310583281642480E+01 0.299668864388672E+01 0.617551331047668E+00 
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TABLE II: The probabilities (in %) of observing selected electron final states in the helium ion 
^He+ and tritium atom ^H arising in the exothermic nuclear reaction Eq.(l) of the three-electron 
lithium-6 atom in its ground ^5"— state with slow neutrons. 



atom/state Is 2s 2p 

'^He+ 9.52386-10-2 1.66605-10"2 6.60170-10-'' 

^H 2.49810-10-3 3.05880-10-'' 2.00251-10-^ 
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